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ABSTRACT 

We present a detailed analysis of the dayside atmosphere of the hot-Neptune GJ 436b, based on 
recent Spitzer observations. We report statistical constraints on the thermal and chemical properties 
of the planetary atmosphere, study correlations between the various molecular species, and discuss 
scenarios of equilibrium and non-equilibrium chemistry in GJ 436b. We model the atmosphere with 
a one-dimensional line-by-line radiative transfer code with parameterized molecular abundances and 
temperature structure. We explore the model parameter space with 10 6 models, using a Markov chain 
Monte Carlo scheme. Our results encompass previous findings, indicating a paucity of methane, an 
overabundance of CO and CO2, and a slight underabundance of H2O, as compared to equilibrium 
chemistry with solar metallicity. The concentrations of the species are highly correlated. Our best-fit, 
and most plausible, constraints require a CH4 mixing ratio of 10~ 7 fol0~ 6 , with CO > 10 -3 , CO2 
~ 10~ 6 iolCP 4 , and H2O < 1CU 4 ; higher CH4 would require much higher CO and CO2. Based on 
calculations of equilibrium and non-equilibrium chemistry, we find that the observed abundances can 
potentially be explained by a combination of high metallicity (~ 10 x solar) and vertical mixing with 
K zz ~ 10 6 — 10 7 cm 2 /s. The inferred metallicity is enhanced over that of the host star which is known 
to be consistent with solar metallicity. Our constraints rule out a dayside thermal inversion in G J 436b. 
We emphasize that the constraints reported in this work depend crucially on the observations in the 
two Spitzer channels at 3.6 Lira and 4.5 Lira. Future observations with warm Spitzer and with the 
James Webb Space Telescope will be extremely important to improve upon the present constraints on 
the abundances of carbon species in the dayside atmosphere of GJ 436b. 

Subject headings: planetary systems — planets and satellites: general — planets and satellites: indi- 
vidual (GJ 436b) — radiative transfer 



1. INTRODUCTION 

The last decade in exoplanetary science has demon- 
strated our capability in detecting and characterizing at- 
mospheres of transiting extrasolar giant planets. Several 
observations have been reported using the Hubble Space 
Telescope (HST), the Spitzer Space Telescope (Spitzer), 
and from ground. Beginning with the first detection 
of sodium in the atmosphere of HD 209458b in tran- 
sit (Charbonneau et al. 2002), and the first detections 
of dayside thermal emission from hot Jupiters TrES-1 
and HD 209458b (Charbonneau et al. 2005; Deming ct 
al. 2005), atmospheric observations of giant exoplanets 
today are a norm, albeit still very challenging (e.g. Knut- 
son et al. 2008; Charbonneau et al. 2008; Swain et al. 
2008; Grillmair et al. 2008; Desert et al. 2009; Swain 
et al. 2009a) . The intensity of observational efforts have 
been matched with equally challenging accomplishments 
in theoretical modeling of cxoplanct atmospheres, and 
data interpretation (Seager & Sasselov, 2000; Seager et 
al. 2005; Burrows et al. 2006 & 2008; Fortney et al. 
2006; Barman et al. 2005; Tinetti ct al. 2007; Show- 
man et al. 2009; Madhusudhan & Seager, 2009). Sev- 
eral inferences of gaseous H 2 0, CH 4 , CO and CO2, and 
thermal inversions have subsequently been made in hot 
Jupiter atmospheres (Burrows et al. 2007; Tinetti et al. 
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2007; Barman, 2007; Grillmair et al. 2008; Swain et al. 
2008; Madhusudhan & Seager, 2009; Swain et al. 2009a 
& 2009b). 

A new era in exoplanetary science has now dawned. 
Latest observations are leading to discovery and char- 
acterization of transiting exoplanets much less massive 
than the archetypal hot Jupiters, namely extrasolar Nep- 
tunes and super-Earths. Several low-mass transiting 
planets are presently known, e.g. hot Neptunes GJ 436b 
(Butler et al. 2004; Maness et al. 2007), HAT-P-llb 
(Bakos ct al. 2010), HAT-P-26b (Hartman et al. 2010), 
and Kepler-4b (Borucki et al. 2010), and super-Earths 
CoRoT-7b (Leger et al. 2009) and GJ 1214b (Charbon- 
neau et al. 2009). These planets have opened a new 
regime in atmospheric modeling and data interpretation. 
However, observations of thermal emission have been re- 
ported for only one of these planets so far, GJ 436b (Dem- 
ing et al. 2007, Stevenson ct al. 2010). 

The first transiting hot Neptune known, GJ 436b, with 
a mass of 22.6 M ffi and radius of 4.2 R ffi , orbits an M 
Dwarf at an orbital separation of 0.03 AU (Butler et al. 
2004; Gillon et al. 2007; Maness et al. 2007). The host 
star has an effective temperature of about 3500 K and a 
metallicity consistent with solar (Torres et al. 2008). 
The average density of the planet is 1.5 g/cc (Torres 
et al. 2008), i.e., similar to the bulk density of Nep- 
tune, an ice giant. At an equilibrium temperature of ~ 
700 K, assuming zero albedo and efficient energy circu- 
lation, the density originally hinted of a hot-ice interior 
(Gillon et al. 2007). More detailed studies of the possible 
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bulk composition of GJ 436b indicate that an additional 
layer of H/He atmosphere would be needed to account 
for the observed radius (Figucira et al. 2009; Rogers & 
Seager. 2009). 

The atmosphere of GJ 436b has been a subject of sub- 
stantial interest in the recent past. Deming et al. (2007) 
and Demory et al. (2007) reported independent detec- 
tions of thermal emission from the dayside of GJ 436b 
in the Spitzer 8 /im IRAC channel. Although no mean- 
ingful inferences about the molecular compositions can 
be drawn from a single data-point, model fits to the 8 
/im flux contrast have favored the interpretation of inef- 
ficient day-night energy redistribution for GJ 436b (De- 
mory et al. 2007; Spiegel et al. 2010, unless there are ad- 
ditional unknown optical absorbers in the atmosphere). 
The models used in Demory et al. 2007 and Spiegel et al. 
2010 had assumed equilibrium chemistry. Given the low 
temperatures of GJ 436b compared to hot Jupiters, equi- 
librium chemistry suggests that the planet's atmosphere 
must be abundant in methane and water vapor, and 
be scarce in carbon monoxide (e.g Burrows and Sharp, 
1999). 

Recent observations have suggested distinct depar- 
tures from predictions of equilibrium chemistry models. 
Pont et al. (2009) reported a transmission spectrum of 
GJ 436b obtained in the 1.1-1.9 /im bandpass using the 
HST NICMOS instrument, but found no significant fea- 
ture in the 1.4 /im water band; the spectrum was flat 
at the 2-tr uncertainties. More recently, Stevenson et al. 
(2010) reported planet-star flux contrasts of the dayside 
atmosphere of GJ 436b in six channels of Spitzer broad- 
band photometry and inferred a deficiency of methane 
in the atmosphere of GJ 436b, using models based on 
Madhusudhan & Seager (2009) and the present work. 
The high planet-star flux contrast observed in the 3.6 
/im IRAC channel was central to the low methane re- 
quirement. Even though our models do not assume equi- 
librium chemistry, the inferred methane mixing ratio of 
~ 10 -7 in a hydrogen dominated atmosphere at ~ 700-fiT 
signals a surprising new regime in atmospheric chemistry 
of extrasolar planets. 

In this work, we report detailed statistical constraints 
on the atmospheric properties of GJ 436b, and explore 
channels of equilibrium and non-equilibrium chemistry 
that might explain the observed chemical abundances. 
We first estimate the atmospheric chemical composition 
and temperature structure at different levels of fit to the 
data, using a 1-D line-by-line radiative transfer model for 
exoplanet atmospheres (Madhusudhan & Seager, 2009). 
Our constraints result from exploring the model param- 
eter space with ~ 10 6 models, optimized using a Markov 
chain Monte Carlo optimization scheme. We then use 
detailed calculations of equilibrium and non-equilibrium 
chemistry, and with different mctallicities, to explain the 
observed constraints on the various molecular species. 

Our results indicate that a high metallicity and non- 
equilibrium chemistry are required to explain the molec- 
ular abundances constrained by the observations. The 
requirement of non-equilibrium chemistry is consistent 
with our findings in Stevenson et al. (2010). The obser- 
vations require the presence of CO and CO2 and a sub- 
stantial depletion of CH4 in the dayside atmosphere of 
GJ 436b. In this work, we find that the constraints on the 
CO and CO2 mixing ratios can be explained by vertical 
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Fig. 1. — Molecular absorption features in Spitzer photometric 
bandpasses. The red dotted lines at the top show the six Spitzer 
bandpasses. The black lines show the extent of absorption fea- 
tures due to the corresponding molecules. The gray curve shows 
a hypothetical model spectrum of GJ 436b based on equilibrium 
chemistry, and the black filled circles show the corresponding in- 
tegrated points in the Spitzer channels. The red filled circles with 
error bars show the observations of GJ 436b reported by Stevenson 
et al. (2010). 

mixing in the atmosphere (with K zz = 10 6 — 10 7 cm 2 /s) 
and a high metallicity (lOx solar). And, we suggest that 
the low CH4 required can potentially be explained by 
non-equilibrium chemistry. We also find that the obser- 
vations require inefficient day-night energy redistribution 
in GJ 436b, which can be confirmed by future observa- 
tions of phase curves. 

In what follows, we first give a brief overview of model 
independent interpretation of Spitzer observations in sec- 
tion^ We then describe our atmosphere model, the pa- 
rameter space optimization scheme, and the models for 
equilibrium and non-equilibrium chemistry we use in this 
work, in section [3] We present the results of our analysis 
in section 2J Finally, in section [SJ we present a summary 
of our work and discuss consequences for future observa- 
tions and theoretical models. 

2. OVERVIEW OF SPITZER DATA INTERPRETATION 

The concentrations of major molecular species in an 
exoplanctary atmosphere can be constrained to some ex- 
tent from observations in the six channels of Spitzer pho- 
tometry, at 3.6 /im, 4.5 /im, 5.8 /im, 8 /im, 16 /im, and 24 
/xm. Figured] illustrates the spectral features of the dom- 
inant molecules in the Spitzer channels, for a synthetic 
thermal spectrum of GJ 436b. CH4 has strong features 
almost exclusively in the 3.6 /im and 8 /im channels. CO 
has a strong feature in the 4.5 /xm channel, also contribut- 
ing to the 5.8 /im channel. And, CO2 has strong features 
in the 4.5 fjm and 16 /im channels. Although H2O has 
several spectral features in most of the Spitzer channels, 
the stronger and/or unique features lie in the 3.6 /im, 5.8 
/im and 24 /im channels. Given the strong features of the 
molecules in one or more Spitzer channels, reasonable 
constraints can be placed with high S/N observations. 
In principle, however, the presence of additional species 
in significant amounts, due to non-equilibrium chemistry 
(e.g. Zahnlc et al. 2009b), and non-LTE effects (Swain 
et al. 2010) can also affect the emergent spectrum. We 
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discuss these factors in section 13.11 

The constraints on the atmospheric properties, that 
might be possible from high S/N Spitzer photometry, can 
be understood in the light of the key parameters effecting 
the emergent spectrum (Madhusudhan & Seager, 2009). 
Under the assumption that H 2 0, CH4, CO and CO2, 
are the dominant spectroscopically active molecules in 
the Spitzer bands, the molecular concentrations consti- 
tute four free parameters. And, although the temper- 
ature structure can involve many free parameters, the 
temperature at the base of the atmosphere (at P ^ 1 
- 10 bar), and the thermal gradient, are the two most 
important parameters, in the absence of a thermal inver- 
sion. Apart from the data themselves, an additional con- 
straint on the parameters appear in the form of energy 
balance (discussed in section l3~Tj) . Thus, the six obser- 
vations can potentially lead to meaningful constraints on 
the six key atmospheric parameters, mainly on the four 
dominant molecules. These same molecular species can 
then provide constraints on the C/H and O/H ratios in 
the atmosphere. 

As an example of constraining molecular abundances 
from the Spitzer data, consider a planetary atmosphere 
with T oq ~ 700 — 1000 K, and a temperature profile de- 
creasing outward (i.e., no thermal inversion). Further 
suppose that a low thermal flux is observed in the 3.6 
/jm channel. Such a low flux would suggest strong ab- 
sorption due to a high concentration of methane. And, 
if the low 3.6 /xm flux is indeed due to a high methane 
abundance (since water vapor also absorbs weakly in the 
3.6 /im channel), the flux should also be low in the 8 //m 
channel where methane also absorbs strongly (as shown, 
for example, by the black filled circles in Figure [1} . On 
the other hand, high planet fluxes in the 3.6 and 8 /im 
channels indicate low absorption due to methane, and 
hence a paucity of methane in the planet atmosphere. 

A second example concerns constraining the molecu- 
lar abundances of CO and CO2. A low flux observed in 
the 4.5 /im channel must indicate atmospheric absorption 
due to CO and/or C02- At 4.5 fxm, CO2 has stronger 
absorption cross-section than CO, so that a small concen- 
tration of CO2 can produce an absorption feature com- 
parable to that from a relatively large concentration of 
CO. The degeneracy between CO and CO2 contributions 
in a 4.5 /jm measurement can be broken by an observa- 
tion in the 16 /im channel, where only CO2 contributes, 
among the two. As a third example, constraints on H2O 
are based primarily on fluxes in the 5.8 /im and 24 /jm 
channels. An important note concerns the presence of a 
thermal inversion in the atmosphere, a region where tem- 
perature increases outward. In such a case, the molecular 
features would be emission features instead of absorption 
features, thus reversing the logic of inferences described 
above (Madhusudhan & Seager, 2010). 

The recent observations of the dayside atmosphere 
of GJ 436b by Stevenson et al. (2010) represent a 
quintessential example of the above inferences. Their ob- 
servations indicate an extremely high flux in the 3.6 /im 
channel and an extremely low flux in the 4.5 /an channel, 
causing a brightness temperature differential of ~ 450 K 
between the two adjacent channels. Thus, following the 
arguments described above, the observations point to- 
wards an extremely low methane abundance and high CO 
and/or CO2 abundances, for a temperature profile with- 



out a thermal inversion. This identification of low CH4 
and high CO and/or CO2 was in fact the central result of 
Stevenson et al. (2010), based on a more elaborate atmo- 
spheric modeling procedure, also discussed in the present 
work. We emphasize that, although the six channels of 
photometry can yield statistical constraints on the atmo- 
spheric properties, they cannot yield a unique solution, 
given the large number of free parameters (Madhusudhan 
& Seager, 2009). 

3. MODEL 

Our goal is to determine the best fitting interpretation 
for observations of the dayside atmosphere of GJ 436b. 
We first fit the data with a large ensemble of 1-D dayside 
atmosphere models of GJ 436b, and determine regions of 
the parameter space that fit the data best. Our results 
yield best-fit constraints on the molecular abundances 
and temperature structure. We then use some of the 
best-fit temperature profiles, along with independent cal- 
culations of equilibrium and non-equilibrium chemistry 
to see if we can explain the observed constraints on the 
molecular abundances. 

3.1. Radiative Transfer Model 

In order to fit the observations with model spectra, we 
use the 1-D exoplanetary atmosphere model developed 
in Madhusudhan & Seager (2009). Our model consists 
of a line-by-line radiative transfer code, with constraints 
of hydrostatic equilibrium and global energy balance, 
and coupled to a parametric pressure-temperature (P-T) 
structure and parametric molecular abundances. This 
modeling approach allows one to compute large ensem- 
bles of models, and to efficiently explore the parameter 
space of molecular abundances and temperature struc- 
ture. 

The major difference of our model from traditional 1-D 
atmosphere models is in the treatment of energy balance. 
Our model requires energy balance at the top of the at- 
mosphere, instead of an iterative scheme to ensure layer- 
by-layer radiative (or radiative + convective) equilibrium 
which is assumed in conventional models (e.g. Seager et 
al. 2005, Burrows et al. 2006, Fortney et al. 2006). For 
a given set of model parameters, we require that the net 
energy output at the top of the atmosphere is less than 
or equal to the net energy input due to the incident stel- 
lar flux; a deficit indicates energy redistributed to the 
night-side. Models where the emergent flux exceeds the 
incident flux are discarded (see Madhusudhan & Seager, 
2009). 

In this work, we consider well mixed atmospheres, i.e. 
uniform mixing ratio of each molecular species over the 
entire atmosphere. The present approach allows us to 
sample a wider range of compositions independent of 
any assumptions of equilibrium chemistry. The molec- 
ular species in our models include molecular hydrogen 
(H2), water vapor (H 2 0), carbon monoxide (CO), carbon 
dioxide (CO2), and methane (CH4). We have used an al- 
ternate approach in previous works (e.g. Madhusudhan 
& Seager, 2009), where we parameterized the abundances 
in terms of deviations from equilibrium chemistry. We 
find that the constraints on the overall mixing ratios do 
not depend critically on the choice of paramctrization. 
Our H2O, CH4, and CO molecular line data are from 
Freedman et al. 2008, and references therein. Our CO2 
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data are from Freedman (personal communication, 2009) 
and Rothman et al. (2005). And, we obtain the H2-H2 
collision- induced opacities from Borysow ct al. (1997), 
and Borysow (2002). 

Our model does not include effects of non-LTE ra- 
diative transfer. Models of exoplanctary spectra with 
Non-LTE radiative transfer calculations have not been 
reported. Recently, Swain et al. (2010) reported detec- 
tion of excess emission at 3.25 fjja in a ground-based ther- 
mal spectrum of HD 189733b, which they surmised to be 
due to non-LTE methane emission. However, follow-up 
observations of Mandell et al. (2010) failed to detect the 
feature reported by Swain et al. (2010). Furthermore, 
their estimates of potential contribution from resonant 
florescence, a non-LTE emission mechanism, of methane 
yielded fluxes too low to significantly contribute to the 
emission spectrum, contrary to the observations of Swain 
et al. (2010). Nevertheless, detailed atmospheric mod- 
els of GJ 436b in the future might need to account for 
non-LTE contributions to spectra, as high-resolution ob- 
servations become possible in the future with the James 
Webb Space Telescope. 

We also do not include in our list of molecules higher 
hydrocarbons that might be potential byproducts of non- 
equilibrium chemistry. Non-equilibrium chemistry can 
lead to hydrocarbons like C2H2 and C2H4, depending on 
the temperature and degree of vertical mixing (Zahnlc 
et al. 2009b; Line et al. 2010). However, as will be dis- 
cussed in section 14.2.21 the observations of G J 436b in- 
dicate high temperatures (T > HOOK) at 1-10 bar pres- 
sures, thereby favoring the oxidation of CH4 to CO, over 
conversion of CH4 to higher hydrocarbons (Zahnle et al. 
2009b). Nevertheless, constraining these species observa- 
tionally would still be a worthwhile exercise when higher 
resolution observations become available with the JWST 

3.2. Parameter Space Exploration for Model Fits 

We use the Markov chain Monte Carlo (MCMC) 
method to explore the model parameter space. The 
MCMC method is a Bayesian parameter estimation algo- 
rithm which allows the calculation of posterior probabil- 
ity distributions of the model parameters conditional to 
a given set of observations, and prior probabilities (see 
e.g. Tegmark et al. 2004; Ford, 2005). In this work, 
our goal is not parameter estimation - the number of 
model parameters (N = 10) exceeds the number of avail- 
able observations (-/Vobs = 6) and rendering the problem 
under-constrained . However, the MCMC method al- 
lows an efficient means of exploring the parameter space 
in search of regions which provide the best fits to the ob- 
servations. We, therefore, use the MCMC method with 
a Metropolis-Hastings scheme within the Gibbs sampler, 
for fine sampling of the model parameter space . And, 
we report error contours in the space of the molecular 
compositions and temperature structure. Our statistic of 
choice is £ 2 , defined as x 2 /N bs (Madhusudhan & Seager, 
2009). In this metric, a £ 2 = 1 indicates models fitting 
the observations within the l-cr observational uncertain- 
ties, on average. Similarly, £ 2 of 2 and 3 indicate fits at 
the 1.41-cr (i.e \^2) and 1.73 -a error bars, respectively. 

Our model described in § 13.11 above has ten free pa- 
rameters (Madhusudhan & Seager, 2009). Six parame- 
ters concern the P-T profile: To, Pi, P2, P3, ai, and a 2 ■ 




500 1000 1500 2000 
T (K) 

Fig. 2.— Pressure-Temperature (P-T) profiles of GJ 436b.The 
purple, red, orange, green and blue profiles correspond to models 
which fit the observations to within § 2 of 1.0, 2.0, 3.0, 4.0 and 5.0, 
respectively. 

And, four parameters correspond to the uniform molecu- 
lar mixing ratios: /h 2 o, /co, /ch 4) and f C o 2 - We define 
the mixing ratio of a molecule as the number fraction 
with respect to molecular hydrogen. 

We define some physically motivated boundaries in the 
parameter space explored by the Markov chain. We 
restrict all the molecular mixing ratios to a range of 
[10 -10 , 0.1]. We also impose the constraint of global 
energy balance by restricting 7/ to [0.0,1.0], where, rj = 
(1 — A)(l — f r ) is the ratio of emergent flux output on 
the day side to incident stellar flux input on the dayside, 
weighted appropriately (Madhusudhan & Seager, 2009). 
Here, A is the Bond Albedo and f r is the day-night en- 
ergy redistribution. The "fit" parameters for the MCMC 
are T , log(Pi), log(P 2 ), log(P 3 ), ai, a 2 , log(/ H2 o), 
log(/co), log(/ C H 4 ), and log(/ C o 2 )- We consider uni- 
form priors in all the parameters. For each system under 
consideration, we run one chain of 10 6 links, which takes 
~ 22 hours on a single processor. 

3.3. Chemistry Model 

After we obtain the constraints on the molecular abun- 
dances from model fits to data, we investigate processes 
of atmospheric chemistry that could explain the required 
abundances. Atmospheres of giant planets in the solar 
system and those of brown dwarfs have revealed the inter- 
play between equilibrium and non-equilibrium chemical 
processes in hydrogen-rich atmospheres (Prinn & Bashay, 
1977; Fegley & Lodders, 1994; Noll et al., 1997; Saumon 
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Fig. 3. — Constraints on the atmospheric properties of GJ 436b. Mixing ratios are shown as ratios by number density. The contours 
show surfaces of minimum £ 2 in the space of atmospheric composition and temperature structure. The purple, red, orange, green and blue, 
correspond to minimum § 2 of 1.0, 2.0, 3.0, 4.0 and 5.0, respectively. The gray and black regions correspond to models which have methane 
mixing ratios greater than 10 — 7 , along with different conditions on the compositions of the remaining molecules, and allowing fits within 
£2 = 3. For the gray surfaces, C0 2 < 10" 3 , H 2 > 10" 5 , and CO < 10" 1 . C0 2 of ~ 10" 3 and CO of - lO" 1 are implausible either in 
equilibrium or non-equilibrium chemistry (Lodders & Feglcy, 2002; Zahnle et al. 2009a), however, we show these solutions for completeness. 
The black surfaces show solutions within more plausible limits (sec section [4TT7TJ of C0 2 < 10~ 4 , H 2 > 10 — 4 , and CO< 10 — 2 . The black 
contour in the C/H - O/H plane requires C/H and O/H > 10 X solar abundances. 
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Fig. 4. — Observations and model spectra for dayside emission from GJ 436b. The black filled circles with error bars show the Spitzer 
observations in the six photometric channels, from Stevenson et al. 2010. The blue dashed lines shows a planet blackbody spectra at 800 
K. The green and red curves show two model spectra , and t he colored circles show the corresponding channel integrated model points. The 
green model is a best-fit model spectrum (see section T4.1. II with non-equilibrium molecular mixing ratios of H2O = 10 — 4 , CO = 7 X 10 -3 , 
CH4 = 10 — 6 , and CO2 = 6 X 10 -6 . The red model has a composition close to chemical equilibrium with solar abundances, with H2O 
= 2 X 10 — 3 , CO = 10~ 5 , and CH4 = 7 X 10~ 4 ; it also contains CO2 = 10 -6 . Both models have the same pressure-temperature profile, 
shown in the inset. The green, non-equilibrium, model has a maximum day-night energy redistribution fraction (/ r ) of 0.03, i.e for zero 
bond albedo (j4g). On the other hand, the red, equilibrium, model has very efficient redistribution, f r < 0.63, for As = 0. 



et al. 2006). At high pressures deep in a planetary at- 
mosphere, molecules react fast enough that all species 
are in thermochemical equilibrium. As the pressures de- 
crease with increasing altitude, thermochemical reaction 
rates decrease, allowing for competing non-equilibrium 
processes with shorter timescales to shift the involved 
species out of equilibrium. We compute the atmospheric 
compositions in equilibrium using the equilibrium chem- 
istry code adapted from Seager et al. (2005). 

3.3.1. Equilibrium chemistry 

Our objective here is to determine whether the ob- 
served constraints on the molecular mixing ratios are 
consistent with chemical equilibrium. To this end, we 
calculate the equilibrium compositions of the species us- 
ing the equilibrium chemistry code originally developed 
in Seager et al. (2000), and subsequently used in Seager 
et al. (2005) and Miller-Ricci et al. (2009). We calculate 
the gas phase molecular mixing ratios for 172 molecules, 
resulting from abundances of 23 atomic species, by min- 
imizing the net Gibbs free energy of the system. The 
multi-dimensional Newton- Raphson method described in 
White et al. (1958) was used for the minimization. 
We adopt polynomial fits for the free energies of the 
molecules, based on Sharp & Huebner (1990). We as- 
sume a hydrogen dominated atmosphere for GJ 436b, 



and we compute equilibrium concentrations of all the 
species at varying metallicites (see section 14. ip , over a 
grid in pressure and temperature. 

At the temperatures of GJ 436b, the most abundant 
and spectroscopically active molecules in the Spitzer 
bandpasses are expected to be H 2 0, CH 4 , CO and CO2. 
Ammonia (NH3) should also be abundant, but does not 
contain major features in the Spitzer channels. At high 
temperatures, e.g. > 1300 K at 1 bar pressure, like those 
in hot Jupiter atmospheres, CO is predicted to be the 
dominant carbon and oxygen bearing species. At lower 
temperatures, on the other hand, CH 4 is the dominant 
carbon bearing species in equilibrium. Water vapor is 
a major carrier of oxygen in either regime. The specific 
amounts of each of these species also depend strongly 
on the assumed metallicity, and pressure. Finally, the 
amount of carbon dioxide in equilibrium is a very strong 
function of metallicity. At the temperatures of GJ 436b, 
solar metallicity yields a CO2 mixing ratio up to 10~ 7 , 
whereas amounts as high as 10~ 4 can be obtained for 30 
x solar metallicity. 

3.3.2. N on- equilibrium chemistry 

Vertical mixing can drive species out of equilibrium in 
regions of the atmosphere where the timescale of vertical 
transport is shorter compared to the timescale govern- 
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ing chemical equilibrium between the relevant species. 
At high pressures, deep in the atmosphere, convection 
is a natural mixing mechanism. However, above the 
radiative-convective boundary, atmospheric instabilities 
and turbulent processes, such as wave breaking, can also 
lead to vertical motions. This form of mixing in the ra- 
diative zone is collectively termed as 'eddy mixing' or 
'eddy diffusion', and is parameterized as a diffusion pro- 
cess, with a coefficient of eddy diffusion (K zz ). Such mix- 
ing shifts molecular species, in radiative regions, away 
from their equilibrium concentrations. Eddy mixing has 
been known in the context of atmospheres of solar system 
giant planets (Prinn & Bashay, 1977; Fegley and Lod- 
ders, 1994), and brown dwarfs (Noll et al. 1997; Saumon 
et al. 2003 & 2006; Hubeny & Burrows, 2007). And, 
Cooper & Showman (2006) studied eddy mixing in hot 
Jupiter atmospheres. Eddy mixing offers a viable ex- 
planation to the excess carbon monoxide discovered in 
hydrogen dominated atmospheres (see e.g Yung & De- 
more, 1999). The primary reaction governing the rela- 
tive abundances of CH 4 and CO in equilibrium is given 
by: 

CO + 3H 2 ^ CH 4 + H 2 0. (1) 

This reaction favors CO at high temperatures and CH4 
at the low temperatures of GJ 436b. Thus, based on equi- 
librium chemistry, it is expected for CO to be dominant 
in hotter lower regions of the atmosphere (T > 1200X), 
and CH4 to be dominant at higher altitudes where tem- 
peratures are lower (T < 1000-fsT). However, atmospheric 
spectra of solar system planets and cooler brown dwarfs 
suggest significant amounts of CO in the upper layers 
of the atmosphere (see e.g. Fegley and Lodders, 1994; 
Noll et al. 1997, Stephens et al. 2009). This is achieved 
by eddy mixing which vertically transport CO from the 
lower regions of an atmosphere to the upper regions. 

Eddy mixing dominates when the mixing time scale 
(Tinix) is shorter than the chemical time scale (r c h cm ) of 
CO in the forward reaction in ([TJ. The forward reaction 
in ([T]) in fact proceeds via multiple steps, and Yung et 
al. (1988) suggested the rate determining step in the 
reaction chain to be (but, cf. Visscher et al. 2010): 

H + H 2 CO + M-> CH3O + M (2) 
The life time of CO is given by: 

[CO] _ [CO] , 
chom ~ d[CO]/dt jfc/ppaCO]' [ > 

where, kf is the rate constant for the forward reaction. 
kf is not known directly, but can be estimated from 
the reverse reaction rate constant (k r ), which is known 
from laboratory experiments, and the equilibrium con- 
stant (K eq ) for the reaction (see Griffith & Yelle, 1999, 
for a detailed discussion). In the present study, we use 
the following estimate of kf, based on Line et al. (2010): 

kf = 3.07 X l -12 T -1.2 e (3927/T) (4) 

Then, r c hem can be calculated from ([3]), using equilibrium 
concentrations of the CO, H and H2CO. 

The mixing time (T m i x ) in radiative regions of the at- 
mosphere is determined by the eddy diffusion coefficient 



(K zz ) and a characteristic length scale for mixing (L), as 

T mix ~ -T? (5) 

L is typically chosen to be the scale height (H) (e.g. 
Prinn & Bashay, 1977; Line et al. 2010). However, for 
the coolest of giant planets, like Jupiter, L can be as low 
as 0.1H (Smith, 1998). In this work, we choose L = 
H; a lower L docs not significantly alter our results, as 
discussed in section |4". 2. II 

As is evident from the above discussion, both T m ; x and 
Tchcm vary with height in the atmosphere; although, r m i X 
varies to a lesser extent than r c i iem . But, while r c i lcm in- 
creases towards higher levels in the atmosphere (i.e with 
decreasing pressure), r m i x increases in the opposite direc- 
tion. The pressure (po) at which r m j X = T c hom is called 
the "quench" level. Above this pressure, i.e in deeper lay- 
ers of the atmosphere, the species are in chemical equi- 
librium, and below this pressure the concentration of the 
species (CO in this case) is fixed, or "quenched" , at the 
equilibrium value at pq. This yields a uniform mixing 
ratio profile for the species for pressures below pq. 

Thus, a higher po implies that CO can be dredged up 
from deeper levels in the atmosphere. And, since in equi- 
librium CO concentration increases with pressure, a high 
Pq implies higher concentration of CO in the upper layers 
of the atmosphere. It can be shown that, for a given r c h e m 
profile, pq increases monotonically with K ZZl thus corre- 
lating a high CO concentration in the upper atmosphere 
with a high K zz (see for e.g. Griffith & Yelle, 1999). 
In the discussion here, we have assumed a temperature 
profile that increases monotonically with pressure. The 
dependence of observed CO on K zz deviates from this 
monotonic behavior for more complicated temperature 
structures, for example in the presence of thermal inver- 
sions or partial isotherms, as will be shown in scction l4~2l 

In this work, we explore the CO mixing ratios result- 
ing from different combinations of K zz and metallicities. 
We explore values of K zz between 10 2 - 10 10 , and metal- 
licities of solar - 30 x solar. Our choices of metallicities 
are motivated by the constraints on CO and CO2 abun- 
dances which indicate high C/H and O/H ratios apriori. 
And, our range in K zz encompasses values found in solar 
system planets and brown dwarfs. For comparison, plan- 
etary atmospheres in the solar system have K zz ranging 
between 10 5 — 10 9 cm 2 /s, and K zz for brown dwarf at- 
mospheres can be as low as 10 2 — 10 4 cm 2 /s (Prinn and 
Bashay, 1977; Yung and Demore, 1999; Saumon et al. 
2003). 

4. RESULTS 

In this section, we present constraints on the atmo- 
spheric properties of the dayside atmosphere of GJ 436b, 
as placed by the six channel Spitzer photometry. We 
first report constraints on the molecular abundances, 
and discuss the correlations between the various species. 
We then discuss the physical plausibility of the solu- 
tions, and present calculations of equilibrium and non- 
equilibrium chemistry attempting to explain the ob- 
served constraints. Finally, we present constraints on the 
atmospheric temperature structure and day-night energy 
redistribution. Our constraints follow from a detailed 
exploration of the model parameter space with ~ 10 6 
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Fig. 5. — Concentrati ons of CH4, H2O, CO, and CO2, predicted by equilibrium chemistry assuming solar and 30xsolar elemental 
abundances (see section I3.3.1H . Each panel shows contours of molecular mixing ratios (i.e. ratio by number density) with respect to 
molecular hydrogen, in pressure-temperature space. The left panels show mixing ratios for solar metallicity, and the right panels ("30xS") 
have 30 X solar metallicity. The gray, brown, red, orange, yellow, green, cyan, blue and purple contours correspond to mixing ratios greater 
than 10 -2 , 10 — 3 , 10 — 4 , 10 — 5 , 10 — 6 , 10 — 7 , 10 — 8 , 10 — 9 , and 10 °, respectively. The black lines are some best fitting pressure-temperature 
profiles. 

by (Madhusudhan & Seagcr, 2009): 

fi,model 



models, using the procedures described in section [ 

Our modeling approach allows placing statistically ro- 
bust model constraints from the data. The pressure- 
temperature (P-T) profiles explored by our models are 
shown in Figure [2] color coded by their degree of fit to 
data. The goodncss-of-fit contours in the space of atmo- 
spheric composition arc shown in Figurc|31 As mentioned 
in section ROI our statistic for the goodness-of-fit is given 
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(6) 



where, fi t model is the planet-star flux contrast of the 
model in each channel, and fi <0 b s and <7j j0 b s are the ob- 
served flux contrast and the 1-a uncertainty, in that 
channel. N b s is the number of observations. Here, 
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N b s = 6, corresponding to the six channels of Spitzer 
photometry. In this metric, a £ 2 = 1 indicates models 
fitting the observations within the 1-a observational un- 
ccrtainitics, on average. Similarly, £ 2 of 2 and 3 indicate 
fits at the 1.41-er (i.e a/2) and 1.73 -a error bars, respec- 
tively. 

The constraints depend on the £ 2 surface one chooses 
for interpretation, apart from any conditions of physi- 
cal plausibility one would like to impose on the models. 
We first report constraints at the £ 2 — 1 and £ 2 = 2 
levels, and with only the barest assumptions of physical 
plausibility. We then discuss additional constraints that 
result from considering some nominal conditions of phys- 
ical plausibility, from equilibrium and non-equilibrium 
chemistry. 

4.1. Constraints on Chemical Composition 

The constraints on the molecular abundances are 
strongly influenced by the correlations between them. 
The correlations between the molecules result from their 
overlapping absorption features in the Spitzer channels 
as described in section [5J We present constraints on the 
mixing ratios of methane (CH4), water vapor (H2O), car- 
bon monoxide (CO), and carbon dioxide (CO^)- Because 
the abundances of molecules arc correlated, constraints 
on any molecule have to be discussed with respect to 
abundances of one or more of the remaining molecules. 

The constraints on all the molecules and the correla- 
tions between them are shown in Figure [3j 

Methane CH4: Our results indicate a substantial 
paucity of methane in the dayside atmosphere of GJ 
436b. Our results place an absolute upper-limit on the 
mixing ratio of methane to be 3 x 1CP 6 — 6 x 10~ 6 , for 
£ 2 ranging between 1-5, and assuming nothing about 
the remaining molecules. However, these upper-limits 
allow for a wide range of abundances of the remaining 
molecules, including some manifestly impractical values. 
Primarily, the constraints include CO2 abundances as 
high 0.3, implying 30% of a hydrogen dominated atmo- 
sphere to be composed of CO2! Assuming a high metal- 
licity for the planet atmosphere (about 30 x solar), CO2 
mixing ratios as high as ~ 10~ 4 can be attained by equi- 
librium chemistry, as shown in Figure [5] (also sec Zahnlc 
et al. 2009a, 2009b). 

The methane mixing ratio is constrained to values be- 
low 10~ 6 , if we impose a plausible limits on the CO2 
abundance. A generous upper-limit on the CO2 abun- 
dance can be assumed to be ~ 10~ 3 , based on the argu- 
ments above. Allowing a maximum CO2 of 10~ 3 , at the 
£ 2 < 1 surface (purple surfaces in Figure [3j) the methane 
mixing ratio is constrained to between 10~ 7 — 10~ 6 , 
for CO2 mixing ratios between 10~ 7 — 10~ 3 . And, at 
the £ 2 < 2 surface (red surfaces in Figure [3]) , CH4 = 
10~ 7 - 10~ 6 , for C0 2 = 10~ 8 - 10~ 3 . There is no lower 
bound on the CH4 abundance; mixing ratios below 10~ 9 
do not have discernible features at the resolution of the 
current data. 

The low methane requirement is enforced primarily by 
the hight planet-star flux contrast in the 3.6 /zm Spitzer 
IRAC channel. Inflating the uncertainties in the 3.6 /an 
channel does not obviate the low methane requirement. 
And, the strong correlation of methane with CO2 arises 
from the large flux differential between the 3.6 /jm and 



4.5 /Ltm channels, as has been described in section [5J 
Methane is also correlated with water vapor which also 
has features in the 3.6 /im channel, albeit to a lesser ex- 
tent, as shown in Figure [3] 

The low mixing ratio of methane is a clear indication 
of non-equilibrium chemistry in the dayside atmosphere 
of GJ 436b, as has been suggested in Stevenson et al. 
(2010). Equilibrium chemistry at the temperatures of 
GJ 436b causes methane to be dominant carbon bearing 
molecule. At solar abundances the methane mixing ratio 
in chemical equilibrium, for typical temperature profiles 
of GJ 436b, is predicted to be 7 x 10" 4 , and 2 x 10~ 2 for 
30 x solar abundances, as is evident from Figure [S] 

Water vapor (H 2 0): Our results place an absolute 
upper-limit on the H2O abundance, as shown in Fig- 
ure [3] The H2O mixing ratio is constrained to < 10 ~ 3 
and < 3 x 10~ 3 for £ 2 < 1 and £ 2 < 5, respectively, if 
we make no assumptions of physical plausibility of the 
solutions. As in the case of methane, however, the H2O 
abundance is also correlated with the CO2 abundance. 
If we restrict CO2 to a generous upper-limit of 10 -3 , as 
described for the case of methane, the £ 2 < 1 and £ 2 < 2 
surfaces constrain the H2O abundance to < 3 x 10~ 4 
and < 10~ 3 , respectively. H2O is also correlated with 
CO and CH4. Thus the H2O abundance can be further 
constrained if we assume conditions of physical plausi- 
bility of all the species simultaneously (discussed below 
in section 14.1. 1[) . The correlations of H2 O with all the 
remaining molecules arise from it numerous features in 
all the Spitzer channels, as described in section [2] 

Carbon monoxide (CO): The abundance of CO is 
highly correlated with the abundance of CO2. If no as- 
sumption is made on the CO2 abundance, the observa- 
tions provide no constraint on the CO mixing ratio, as 
shown in FigureOU If the mixing ratio of CO2 is restricted 
to the conservative upper-limit of 10 -3 , the observations 
constrain CO to be > 3 x 10~ 5 for £ 2 < 1, while still 
leaving it unconstrained at £ 2 ~ 2 and higher. 

Despite the constraints on CO above, simultaneous 
conditions of physical plausibility on all the molecules 
indicate that a very high CO abundance (> 10 -3 ) is es- 
sential to explain the observations. This will be discussed 
in section 14.1.11 below. The strong correlation between 
CO and CO2 arises from the fact that both molecules 
have strong features in the 4.5 /im IRAC channel, the 
CO2 feature being stronger, as described in section [2j 
Although not apparent in Figure [3j CO is also correlated 
with H2O and CH4, via their correlations with C02- 

Carbon dioxide (CO2): The observations require a 
clear presence of CO2 in the atmosphere, as shown in Fig- 
ure HI Since CO2 and CO are correlated, a low concen- 
tration of CO2 requires a high concentration of CO. At 
the £ 2 < 1 surface, a CO2 concentration less than 10 ~ 7 
requires a CO concentration greater than 10~ 2 . On the 
other hand, having a CO2 concentration of 10~ 4 allows 
for CO concentrations as low as ~ 10~ 3 , at the £ 2 < 1 
surface. For £ 2 surfaces of 2 and higher, lower CO abun- 
dances can fit the data for a given CO2 abundance, as 
shown in Figure [3] However, CO2 is also correlated with 
H2O, such that an H2O abundance greater than 3 x 10~ 4 
would require a CO2 abundance greater than ~ 10 -4 , at 
the £ 2 < 1 surface. 

The observations themselves do not place any upper- 
limit on the CO2 abundance. However, the maximum 
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amount of CO2 possible can be constrained based on the- 
oretical limits of equilibrium and non-equilibrium chem- 
istry. For the best-fit temperature profiles of GJ 436b 
shown in Figure [2j and pressures in the 10 ~ 3 — 1 bar 
range, thermochemical equilibrium can yield CO2 mix- 
ing ratios up to 10 for solar metallicity and up to 10 
for ~ 30 x solar metallicities (see Figure [5]). The strong 
dependence of CO2 abundance on metallicity has been 
reported before (Lodders & Fegley, 2002; Zahnle et al. 
2009a) 

4.1.1. Plausibility of the Abundance Constraints 

The constraints reported above assume nothing with 
regards to the physical plausibility of the models, except 
the conservative limit of CO2 < 10 -3 used for sake of ar- 
gument. However, reasonable theoretical constraints can 
be placed over the observed constraints based on well es- 
tablished arguments of equilibrium and non-equilibrium 
chemistry (see section |3~3"]) . The mixing ratios of CH4, 
H2O, CO and CO2 under chemical equilibrium are shown 
in Figure [5j for a range of temperatures and pressures 
pertinent to GJ 436b, along with some best-fit P-T pro- 
files. A detailed discussion of non-equilibrium chemistry, 
via eddy mixing, is presented in section 14.21 

We find that stringent constraints on the molecular 
abundances required by the data can be placed even 
with modest assumptions of atmospheric chemistry. As 
alluded to in section 14.11 above, and shown in Figure 03 
CO2 mixing ratios up to 10 -4 are allowed for high metal- 
licity (also see Zahnle et al. 2009a; Lodders & Fegley, 
2002). H2O is a major carrier of oxygen in the desired 
temperature range. The H2O mixing ratio is expected 
to be ~ 10~ 4 and ~ 5 x 10~ 3 , for solar and 30x solar 
abundances, respectively And, while there is no lower 
limit on the CO abundance, the CO upper-limit is fixed 
by the metallicity; CO < 10 -4 and < 10~ 2 , for solar 
and 30 x solar abundances, respectively. Finally, the con- 
centration of methane in equilibrium follows the carbon 
abundance. At the temperatures of GJ 436b, methane is 
supposed to be highly abundant in equilibrium, as shown 
in Figure[SJ with mixing ratios of 7 x 10~ 4 , and 2 x 10~ 2 
for solar and 30 x solar abundances, respectively. How- 
ever, CO can be enhanced and methane can be depleted 
to some extent due to non-equilibrium chemistry (Zahnle 
et al. 2009b), which will be discussed in section l4~2l be- 
low. 

The constraints due to the considerations of physi- 
cal plausibility are shown in Figure [3J The gray sur- 
face shows regions assuming conservative boundaries of 
C 2 < 3, CH 4 > 10~ 7 , C0 2 < 10" 3 , H 2 > 10^ 5 , and 
CO< 10 _1 . And, the black surfaces show a subset of the 
gray surface with £ 2 < 3, CH 4 > 10~ 7 , C0 2 < 10~ 4 , 
H 2 > 10~ 4 , and CO < 10~ 2 . The black contours rep- 
resent our most likely interpretation of the observations, 
which will be justified below. A best-fit model consistent 
with our constraints above is shown in Figure HJ 

Our most plausible constraints on the atmosphere of 
GJ 436b indicate the possibility of high metallicity, along 
with non-equilibrium chemistry. The black surfaces in 
Figure [3] show that a CO abundance > 10~ 3 is required 
to have a H2O abundance of > 10~ 4 and a methane 
abundance above 10 -7 . The corresponding constraint 
on CO2 is 10~ 6 — 10~ 4 . While this CO2 abundance can 
be explained based on equilibrium chemistry with high 



metallicity alone (see Figure [5]), the high CO abundance 
requires a high metallicity along with non-equilibrium 
chemistry. As shown in Figure [5j very high CO abun- 
dances can exist at the bottom of the atmosphere, for 
high metallicity, however the cooler upper layers of the 
atmosphere have much lower CO abundance. Non- 
equilibrium chemistry in the form of eddy mixing can 
transport CO from the lower layers to the upper layers 
of the atmosphere to cause a uniformly high CO over the 
entire atmosphere, as discussed in section [3.3.21 On the 
other hand, the low CH4 abundance could potentially be 
caused by non-equilibrium chemistry as well (Zahnle et 
al. 2009b). In what follows, we will explore the realm of 
non-equilibrium chemistry in an attempt to explain the 
observed constraints. 

4.2. Explanations for Non- equilibrium Abundances 

Our best-fit constraints on the chemical composition of 
GJ 436b require substantial deviations from equilibrium 
chemistry with solar metallicity. Non-equilibrium pro- 
cesses have been known to influence chemical composi- 
tions of several planetary and brown dwarf atmospheres, 
as discussed in section [3751 Here, we explore channels of 
non-equilibrium chemistry in search of potential explana- 
tions to the observed constraints which are inexplicable 
by equilibrium chemistry alone - namely the high CO 
and CO2 abundances, and the low CH4 abundance. 

4.2.1. High CO and C0 2 

The high abundances of CO and CO2 required by the 
observations can be achieved via eddy mixing, along with 
a high metallicity. As explained in section I3.3.2[ eddy 
mixing transports CO from the deeper layers of the at- 
mosphere to the upper layers of the atmosphere where 
CO is less abundant. The relevant quantity is the quench 
pressure (po), which denotes the pressure level in the at- 
mosphere below which the rate of vertical mixing is faster 
than the reaction rate for chemical equilibrium; the CO 
concentration in the upper atmosphere (P < po) freezes 
at its value at P = p . Figure [JJ] shows the dependence of 
the quench level (po) on the diffusion coefficient (K zz ), 
for a best-fit P-T profile of GJ 436b. Higher values of 
K zz lead to mixing of species from deeper levels of the 
atmosphere, i.c higher po. It can be seen from the figure 
that K zz values between 10 6 - 10 7 can cause mixing from 
a quench level of ~ 1 bar. 

The abundance of CO in the upper atmosphere de- 
pends on K zz and the metallicity. The right panel of 
Figure [JJ] shows the range of CO mixing ratios that are 
attainable with different values of K zz and metallicities. 
For solar metallicity, a maximum CO mixing ratio of 
~ 10~ 4 is attainable for K zz values between ~ 10 7 — 10 8 . 
However, CO mixing ratios above 10 -3 require metal- 
licities greater than 10 x solar, and K zz ~ 10 5 — 10 8 . 
For a very high metallicity of 30 x solar, a maximum 
CO mixing ratio 10 -4 can be attained for a rather low 
K zz of ~ 10 3 cm 2 /s, and a maximum of ~ 10 -2 for K zz 
between ~ 10 6 — 10 s cm 2 /s. For comparison, planetary 
atmospheres in the solar system have K zz ranging be- 
tween 10 5 — 10 9 cm 2 /s (see e.g. Prinn and Bashay, 1977; 
Fegley & Lodders, 1994; Yung and Demore, 1999), and 
for brown dwarf atmospheres, K zz can range between 
10 2 - 10 6 cm 2 /s (Saumon et al. 2006; Stephens et al. 
2009). 
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Fig. 6. — Vertical eddy diffusion in GJ 436b. The two panels show the dependence of the "quench" pressure (po) and the CO mixing 
ratio on the coefficient of eddy diffusion, K zz . "S", 10 X S, and 30 X S refer to solar, 10 times solar and 30 times solar metallicities. 
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Fig. 7. — Non-equilibrium chemistry in GJ 436b. Panel A shows the pressure-temperature (P-T) profile of a best-fit model. Panel B 
shows the times scales of eddy mixing (r m i x ) and of CO-CH4 equilibrium chemistry (T c hem)- Panel C shows the influence of eddy mixing 
on the mixing ratios of the four prominent species. The solid lines in panel C show the compositions from equilibrium chemistry, which 
continue as dotted lines above the quench level, and fixed to the values at the quench level. The dashed lines show the mixing ratios from 
Zahnle et al. (2009b), for an isothermal atmosphere with T = 1200-K", K zz = 10 7 cm 2 /s, and 5xsolar metallicity, which yields a CH4 
mixing ratio below 10 -6 at the P = 1 bar level. The y-axis in all the panels is pressure in bars. The horizontal black dashed line in all the 
panels show the quenching pressure level. 



The constraint on K zz is also dependent on the choice 
of the characteristic length scale for mixing (L), through 
equation [S] In the above results, we used the typical 
assumption of L = H, the scale height of the atmosphere 
(e.g. Prinn & Bashay, 1977; Line et al. 2010). However, 
for the coolest of giant planets, like Jupiter, L can be 
- 0.1H (Smith, 1998). We find that a choice of L < 
H for GJ 436b results in only a modest change in the 



required K zz . For instance, the peak CO mixing ratio 
for 10 x solar metallicity, for the L = H case is attained 
for K zz — 10 6 — 10 s cm 2 /s, where as the same value 
for even the extreme case of L = 0.1H is attained for 
K zz = 10 4 — 10 6 cm 2 /s. Such differences in K zz can be 
compensated by uncertainty in the metallicity, as seen in 
Figure [H 

Finally, a high metallicity is also consistent with the 
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high CO2 abundance required to explain the observa- 
tions. At the quench levels of ~ 1 bar, the equilibrium 
mixing ratio of CO2 is 10 -7 and 10 -4 , for solar and 30 
x solar metallicity, as shown in Figure [5j The strong 
dependence of CO2 on metallicity has been reported is 
previous studies (Lodders & Fegley, 2002; Zahnlc ct al. 
2009a), suggesting CO2 as a key metallicity indicator. 
We, thus, hnd that a metallicity of 10 - 30 x solar, and 
vertical mixing of K zz ~ 10 5 — 10 s , can simultaneously 
explain the constraints on CO and CO2 observed in sec- 
tion EXT1 The K zz and metallicity are constrained fur- 
ther in section 14.2.31 

4.2.2. Low CHi 

The low methane abundance required by the obser- 
vations cannot be explained by equilibrium chemistry. 
At P ~ 1 bar and T ~ 1000 K (corresponding to the 
best-fitting P-T profiles), equilibrium chemistry yields 
methane mixing ratios of about 10 -4 and 10~ 2 , for solar 
and 30 x solar metallicity, respectively, as shown in Fig- 
ure[5] Furthermore, considering non-equilibrium thermo- 
chemistry using the CO-CH4 reaction pathway alone, as 
investigated in section 14 . 2 . 1 1 above . does not explain the 
low methane abundances. By non-equilibrium thermo- 
chemistry, we mean departures from equilibrium chem- 
istry due to vertical mixing. The remaining alternatives 
in searching for a low-methanc solution include a more 
detailed treatment of non-equilibrium thermochemistry, 
including all the possible species and reactions, and pho- 
tochemistry. Zahnlc ct al. (2009b) have reported such 
calculations over a range of temperatures and eddy mix- 
ing coefficients that are relevant to the current situation. 

Substantial depletion of methane below equilibrium 
values is possible due to non-equilibrium thermochem- 
istry and photochemistry. Full kinetics models yield a 
substantial amount of free radicals, especially H (Liang 
et al. 2003; Zahnle et al. 2009b). Zahnlc ct al. (2009b) 
show that the resultant overabundance of H leads to de- 
piction of methane via two main channels. At high tem- 
peratures (T > 1200 K), H readily reacts with H2O to 
yield the reactive OH radical, which oxidizes CH4 in very 
short timescales. The result is an excess production of 
CO. At lower temperatures (T < 1200 K), water is more 
stable, leading to a reducing environment with abundant 
free H. Methane, being still reactive, reacts with H to 
form higher hydrocarbons, and can be markedly reduced 
below equilibrium levels. Thus, in either scenario, the 
methane concentration at the observable pressure levels 
can be substantially depleted below equilibrium, depend- 
ing on the temperature, eddy mixing coefficient (K zz ), 
metallicity, and photochemistry (Zahnle et al. 2009b). 

Two key parameters governing non-equilibrium 
methane chemistry are the temperature and K zz . The 
best-fitting P-T profiles required by the observations 
all have temperatures over 1100 K in the lower atmo- 
spheres (P > 10 bar), similar to those obtained from 
self-consistent models of GJ 436b (Spiegel et al., 2010). 
Figure [7] shows a sample best- fit P-T profile (panel A) , 
which has an isothermal temperature structure of 1200 
K below P ~ 1 bar. And, as described in section 14.2.11 
the observational constraint on the CO concentration 
requires K zz ~ 10 7 cm 2 /s. Panel B of Figure [7] shows 
the CO-CH4 reaction timescale and the eddy mixing 
time scale, for K zz of 10 cm 2 /s, varying with pressure; 



the time scales intersect at the quench pressure (po) of 
~ 1 bar. 

The methane mixing ratios attainable via kintctics and 
photochemistry for the required parameters are discussed 
in Zahnle et al. (2009b). The parameters are: T = 1200 
K, K zz = 10 7 cm 2 /s, and p = 1 bar. Panel C of Figure[7] 
shows mixing ratio profiles from our CO-CH4 disequilib- 
rium model (section l4.2.1[) . and the mixing ratio profiles 
from Zahnle et al. (2009b), corresponding to the same 
T and K zz . As shown in Panel C, their results show 
that CH4 < 10~ 7 is possible at the 1 bar level, assum- 
ing 5 x solar metallicity and a stellar irradiation that is 
100 x the solar insolation at Earth. Lower temperatures 
and higher K zz both favor higher methane concentration 
at the 1 bar level. And, the chemistry is less sensitive 
to photochemistry at higher pressures deep in the atmo- 
sphere. 

4.2.3. Joint Constraints on Metallicity and Eddy Mixing 

The observational constraints on all the molecules yield 
a plausible set of constraints on the metallicity and K zz . 
A full kinetic and/or photochemical model is beyond the 
scope of the current work. However, based on the mod- 
els of Zahnle et al. (2009b) with the T and K zz that 
we constrain for GJ 436b, we find that the low methane 
mixing ratio (10~ 7 — 10~ 6 ) observed can most likely be 
explained by non-equilibrium chemistry. Although such 
low methane abundances can be obtained for 5 x solar 
metallicity and K zz ~ 10 7 cm 2 /s, a higher metallicity 
would be required to simultaneously explain the high CO 
and CO2 abundances, as described in section 14.2.11 On 
the other hand, too high of a metallicity, of say 30 x 
solar, might also increase the CH4 abundance to above 
the favorable levels. We therefore choose an intermedi- 
ate value of 10 x solar which satisfies the CO and CO2 
constraints. 

Based on the above reasoning, the observed chemistry 
in the dayside atmosphere of GJ 436b likely results from 
a high metallicity (~ 10 x solar) and a K zz ~ 10 6 — 10' 
cm 2 /s. A more robust conclusion would be possible with 
a full kinetic + photochemical model of GJ 436b, using 
our best-fit P-T profile, and the appropriate UV spec- 
trum, or flux scaling, for the host star. The present con- 
straints indicate a significant enhancement in the metal- 
licity of G J 436b over that of the host star whose metal- 
licity is known to be consistent with solar (Torres et al. 
2008). Additionally, the observations constrain the C/O 
ratio between 0.5 and 1.0, with the most likely solutions 
(black surfaces in Figure [3]) suggesting a C/O between 
^0.85 and 1. However, whether the methane depiction 
required by the data can be obtained in a high C/O en- 
vironment, e.g. of C/O = 1, needs to be investigated by 
a full non-equilibrium chemistry model in the future. 

4.3. Temperature Structure and Day-Night Energy 
Redistribution 

The Spitzer observations provide important constraints 
on the vertical thermal gradient and the energy balance 
in the dayside atmosphere of GJ 436b. The observed 
brightness temperatures in the six Spitzer channels range 
from about 700 K in the 4.5 (im IRAC channel (3-er upper 
limit) to about 1100 K in the 3.6 fim channel (Stevenson 
et al. 2010). With the exception of the 3.6 ^m IRAC 
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channel, observations in all the remaining five channels 
are consistent with a black-body planet spectrum at 750 
K ± 100 K. However, the 3.6 /jm observation of ~ 1100 
K brightness temperature is a major exception, requir- 
ing a temperature differential of ~ 400 K between the 3.6 
fun and 4.5 /im channels, implying a very steep temper- 
ature gradient in the atmosphere. The dayside pressure- 
temperature (P-T) profiles of G J 436b constrained by the 
observations are shown in Figure [2] The best fitting P- 
T profiles (in purple) have temperatures varying by over 
~ 400 K per bar of pressure in the lower atmosphere, 
required primarily by the large temperature differential 
described above. 

The observations rule out the presence of a significant 
thermal inversion in the dayside atmosphere of GJ 436b. 
Our results show that the observations cannot be fit 
with an inversion model for any chemical composition 
(although very small inversions which might cause only 
weak observable features cannot be ruled out by present 
data) . A significant thermal inversion in this atmosphere 
would have caused the brightness temperatures in the 
4.5, 5.8, and 8 /im channels to be markedly higher than 
the 3.6 /im channel, much higher than what the cur- 
rent data indicate. Finally, as has been known for all 
hot Jupitcrs, the data indicate that the observable day- 
side atmosphere of GJ 436b is mostly radiative, with the 
radiative zones of some of the best fitting P-T profiles 
extending to pressures above ~ 10 bar. The isotherms 
at the high pressure ends of the P-T profiles are sugges- 
tive of the radiative diffusion approximation in the high 
optical depth limit (see e.g. Madhusudhan & Seager, 
2009), and are also found in other self-consistent models 
reported in the literature (e.g. Spiegel et al. 2010). 

The large brightness temperature observed in the 3.6 
(im channel, of 1100 K, is also indicative of low day-night 
energy redistribution in GJ 436b. The 3.6 /jm channel 
probes levels deep in the atmosphere (around pressures 
of 1 bar or higher). A high 3.6 (im brightness temper- 
ature, therefore, indicates a high blackbody continuum 
emerging from the base of the dayside atmosphere at ~ 
1 bar. Our best-fit models show that the net emergent 
flux on the dayside nearly balances the incident stellar 
flux, implying that very little energy is circulated to the 
night side. The bottom-right panel of figure [3] shows 
the £ 2 contours in the rj = (1 — A)(l — f r ) and T e ff. 
n = (1 — A)(l — f r ) is obtained from energy balance, 
where A is the bond albedo and f r is fraction of energy 
redistributed to the night side (Madhusudhan & Seager, 

2009) . 

The best fitting models favor r/ > 0.75, or a maximum 
day-night redistribution (f r ) of 0.25, i.e. for A = 0. An 
albedo of, say 0.1, further restricts the distribution to 
0.2. Our results support similar conclusions arrived at 
by previous works (Deming et al. 2007; Spiegel et al. 

2010) . It is to be noted that we do not assume uniform 
illumination of the planetary dayside by the stellar irra- 
diation (i.e weighing the stellar flux by / = |). Instead, 
we use / = f i according to the prescription of Burrows 
et al. (2008), which is also used in Madhusudhan & Sea- 
ger (2009). Consequently, our estimation of T e g for a 
given r\ is typically higher than one would obtain using 
the / = 5 assumption. 

Our conclusion of a low day-night redistribution on 



this planet assumes special significance for potential fu- 
ture observations of thermal phase curves of GJ 436b. 
Thermal phase curves in the 3.6 /jm and 4.5 /im IRAC 
channels, feasible with warm Spitzer, should show clear 
model-independent evidence of a high day-night temper- 
ature contrast, according to our present results. A find- 
ing on the contrary, i.e finding efficient redistribution in 
the phase curves, can imply the possibility of a substan- 
tial interior energy source in GJ 436b. 

5. DISCUSSION AND SUMMARY 

We presented a detailed analysis of the dayside atmo- 
sphere of GJ 436b. Our results show that a high metal- 
licity along with non-equilibrium chemistry are required 
to explain the observations. We also studied the cor- 
relations between the various molecular species, and re- 
ported detailed constraints on the metallicity, chemical 
processes, and day-night energy circulation. Although 
our results come from observations in six channels of 
Spitzer photometry, some channels arc more important 
than others. Here, we discuss the relative importance 
of the different Sptizer channels to our conclusions. We 
also discuss some potential alternate interpretations of 
the data, followed by a summary of our results. 

5.1. Sensitivity of Results to Spitzer Observations 

The constraints reported in this work depend critically 
on the two Spitzer IRAC observations at 3.6 /im and 4.5 
/im. The high planet-star flux contrast in the 3.6 /im 
channel is responsible for the constraints of low methane 
abundance and low energy circulation. The low flux con- 
trast in the 4.5 /jm channel is responsible for the require- 
ment of high CO and/or C0 2 in GJ 436b. While the 
observation in the 3.6 /im channel was reported to be of 
the higher S/N of all channels, the 4.5 /im channel was 
a non-detection (Stevenson et al. 2010). Nevertheless, 
future observations in both these channels would be ex- 
tremely important to confirm that these fluxes actually 
represent the steady state atmosphere in GJ 436b. The 
observations in the remaining four channels (5.8 /im - 
24 /im) are much less constraining, although still very 
useful. For instance, the moderate flux observed in the 
8 /im channel, where methane absorbs strongly, is im- 
portant to the conclusion that the very high flux in the 
3.6 ^m channel could not have been due to a thermal 
inversion causing methane emission. 

Our constraints on the molecular abundances are fairly 
robust with respect to the observational uncertainties. 
The high flux in the 3.6 /im channel cannot be explained 
by equilibrium chemistry. A model with equilibrium 
chemistry and solar or 30 x solar metallicities predicts 
planet-star flux contrasts that are lower than the ob- 
served value by over 4 -a (also in agreement with mod- 
els of Demory et al. 2007; Spiegel et al. 2010). On 
the other hand, the contrasts predicted in the 4.5 /im 
channel based on equilibrium chemistry alone would be 
higher than the observed non-detection by over 3-cr (also 
see Spiegel et al. 20010). We have been conservative 
in our analysis by allowing our best fits to the 4.5 /im 
point to lie within the 3-cr ± 1-er upper-limits. Had we 
considered this point to be a strict non-detection at 1-a, 
our results would predict even higher CO and/or CO2. 
Finally, the large uncertainties in the observed fluxes in 
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the 16 /^m and 24 /im channels provide only fiducial con- 
straints on the temperature structure and the H2O and 
CO2 abundances. 

5.2. Alternate Interpretations 

The high planet-star flux contrast observed in the 3.6 
/im IRAC channel is central to most of the constraints 
reported in this work. Our inferences could partly be 
restricted by the specific choices that are inherent to our 
models. For example, our models do not include clouds 
or hazes. Although scattering from hazes has been sug- 
gested to be potentially relevant in the optical and near- 
IR (e.g. Sing et al. 2009), a high contribution at the 
longer wavelengths of the 3.6 (im channel is unlikely. We 
have also assumed the planet atmosphere to be in local 
thermodynamic equilibrium (LTE). Swain et al. (2010) 
reported a potential signature due to non-LTE methane 
fluorescent emission in HD 189733b observed in the range 
of the 3.6 /im channel (but c.f. Mandell et al. 2010, who 
did not find such a feature in follow-up observations). 
The dayside atmosphere of GJ 436b might also be vari- 
able, as has been reported previously for hot Jupitcrs 
(Grillmair et al. 2008; Madhusudhan & Seagcr, 2009). 
However, for variability to explain the observed flux in 
the 3.6 /im channel, the temperature at the ~ 1 bar level 
in the atmosphere of G J 436b has to exhibit fluctuations 
greater than 400 K, between subsequent observations of 
Stevenson et al. (2010). These and other alternate expla- 
nations are worth exploring with more data at different 
epochs. 

5.3. Summary 

We have presented constraints on the chemical com- 
position and temperature structure of the dayside atmo- 
sphere of hot Neptune GJ 436b, based on recent Spitzer 
observations. One of our key findings is the strict up- 
per limit on the mixing ratio of methane. We find that 
models fitting the observations require a methane mix- 
ing ratio below 10 -6 . Slightly higher methane mixing 
ratios require CO2 ~ 10 -3 — 10~ 2 , which is implausi- 
ble in the hydrogen rich atmosphere with the tempera- 
ture structure of GJ 436b. The abundances of all the 
molecules are highly correlated. Applying nominal con- 
ditions of physical plausibility, we find the constraints on 
the molecular mixing ratios to be CH 4 - 10~ 7 - 10~ 6 , 
CO > IO- 3 , C0 2 ~ 10~ 6 - 10~ 4 , and H 2 < 10~ 4 . 
These constraints on the molecular abundances cannot 
all be explained based on equilibrium chemistry, for any 
metallicity, as reported in Stevenson et al. (2010). At 
the temperatures of G J 436b, equilibrium chemistry with 
solar abundances predicts CH4, CO and CO2 mixing ra- 
tios to be ~ 5 x 10 -4 , ~ 10 -5 and ~ 10 -7 , respectively, 
contrary to the observed abundances. 

The observed constraints on the molecular abundances 
can be explained by a combination of high metallicity 
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and non-equilibrium thermochemistry. A high metallic- 
ity is required for a high CO2 abundance. Vertical mixing 
along with a high metallicity is required to dredge up the 
high CO abundance from the lower layers of the atmo- 
sphere to observable layers. Finally, vertical mixing and 
photochemistry can cause substantial depletion of CH4, 
as reported by Zahnle et al. (2009b). At the tempera- 
tures (T > 1100 K) we obtain for the lower atmosphere 
of GJ 436b, the results of Zahnle et al. (2009b) indicate 
the depletion of equilibrium CH4 via oxidation to CO, 
caused by overabundance of the H radical. Our joint 
analysis of the parameters of non-equilibrium chemistry 
required to explain the abundances of all the species, sug- 
gests that the dayside atmosphere of GJ 436b has a high 
metallicity of ~ 10 x solar and a diffusion coefficient of 
K zz = 10 6 — 10 7 cm 2 /s. The metallicity is substantially 
enhanced over that of the host star which is consistent 
with solar metallicity (Torres et al. 2008). 

Our results also constrain the dayside temperature 
structure and the day-night energy redistribution in the 
atmosphere of GJ 436b. A temperature inversion is 
ruled out by the current observations; although, small 
inversions which arc not observable at the resolution 
of the current photometric data cannot be conclusively 
ruled out. The observations also suggest inefficient day- 
night energy redistribution (f r ) in GJ 436b, requiring 
(1 - A B )(1 - f r ) = 77 > 0.7, at the 1-a fit. Thus, the 
maximum f r allowed by the data at the 1-a fits is 0.3 
for zero bond albedo (Ab), and 0.23 for Ab = 0.1. Fu- 
ture observations of thermal phase curves in the available 
warm Spitzer channels will be instrumental in validating 
the low redistribution requirement. A finding on the con- 
trary might indicate a substantial interior energy source. 
We emphasize that the constraints reported in this work 
depend primarily on the two Spitzer channels (3.6 /im 
and 4.5 /im). Thus, future observations in these channels 
will be extremely important in confirming the present 
and previous results on the atmosphere of GJ 436b. 

The atmosphere of the hot-Neptune GJ 436b presents 
new challenges and opportunities for detailed modeling 
of exoplanet atmospheres. As low-mass transiting plan- 
ets continue to be discovered, unexpected findings are 
likely to continue. The next generation of models and 
observations will help unravel those mysteries, and help 
put our solar system in perspective. 
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